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SUMMARY 

The available information on the artificial radiation belt formed by 
the July 9, 1962, high altitude nuclear explosion is reviewed. Data 
from Injun (1961 o2), Telstar I (1962 del), Traac (1961 0-772), and 
Ariel I (1962 ol) are combined to form one picture of the artificial 
belt. The data are consistent to about a factor of 3. The flux map ob- 
tained in this way is used to calculate the flux encountered by several 
satellites. These show reasonable agreement with data on solar cell 
damage. Preliminary data on particle lifetimes are presented. Parti- 
cles at L > 1.30 are expected to last several years on the basis of 
coulomb scattering. Crude calculations of shielding are made to indi- 
cate the doses received inside various vehicles. 
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THE ARTIFICIAL RADIATION BELT 
MADE ON JULY 9, 1962* 

by 

Wilmot N. Hess 
Goddard Space Flight Center 


INTRODUCTION 


On July 9, 1962, at 0900:09 UT a nuclear explosion of about 1.4 megatons was carried out at 400 
kilometers above Johnston Island in the Pacific Ocean. This explosion produced, as was expected, an 
artificial radiation belt. However, the intensities in this radiation belt are considerably higher than 
were expected. Three days after the explosion the U.S.-U.K. joint satellite Ariel I (1962 ol) stopped 
transmitting. On August 2, Transit IV-B (1961 on) 1) stopped transmitting; Traac (1961 ar\2) stopped 
on August 14. Instruments on Ariel I, Traac, and Injun (1961 o2) showed large particle fluxes shortly 
after the explosion. It took about a month to start getting some grasp of the characteristics of the 
new radiation belt. This is a status report on the new belt as of September 12. 


AVAILABLE DATA 

The information that is available to form a picture of the new radiation belt comes mostly from 
particle detectors on the Ariel I, Injun, Traac, and Telstar I (1962 ael) satellites. In addition to 
these data we can use the observed solar cell damage on satellites as an integral measurement of the 
trapped electron flux. Also, some data are available from dosimetry measurements. 

Some of the original data about the enhanced trapped particle fluxes after the July 9 explosion 
came from the x-ray detector on the Ariel I satellite (private communication from A. Willmore, 
University College, London). This instrument was not designed to count charged particles and there- 
fore its efficiency is uncertain. The data from it are quite useful in studying the time decay of the 
trapped flux and in locating contours of constant flux in B-L space. 

Data received by the shielded 213 GM counter on Injun have been analyzed to give the first picture 
of the new radiation belt (Reference 1). This counter is the background channel of the magnetic spec- 
trometer, SpB. It has 3-1/2 gm/cm 2 of Pb shielding and about 1 gm/cm 2 of wall and miscellaneous 
shielding. It was supposed to give the penetrating background to be subtracted from the other channels 

•This is an abridged version of Technical Memorandum X-788, a confidential report entitled "The Artificial Radiation Belt. This version 
will also be published in the Journal of Geophysical Research. 
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of the spectrometer. This detector is now called on to provide quantitative information, and it has 
been calibrated after the fact. It is nearly omnidirectional. Fluxes are obtained from the count rates 
by dividing by G 0 = 0.11 cm 2 . Other detectors on Injun also give useful data sometimes, but often 
they are saturated and not usable. So far, little data have been analyzed from any Injun detectors ex- 
cept SpB. 

Telstar I has on it a solid state p-n junction detector with pulse height analysis that selects 
electrons in different energy ranges from 0.2 to 1 Mev (private communication from W. Brown, Bell 
Telephone Laboratories). A lot of data have been reduced from Telstar I for two channels of the 
electron detector. This detector has given all the data currently available at high altitudes. It is 
directional, with an aperture half-angle of about 10 degrees. The fluxes are made omnidirectional 
by multiplying by the appropriate solid angle factor and then using a factor between 1 and 2 to cor- 
rect, roughly, for the nonisotropic angular distribution. 

Traac has a 302 GM counter shielded by 0.265 gm/cm 2 of Mg, which will count electrons of 
energy above 1.5 Mev (private communication from G. Pieper and L. Frank, Applied Physics Labora- 
tory). It is essentially omnidirectional. Fluxes are obtained by dividing by G o = 0,75 cm 2 and cor- 
recting for saturation for high count rates. 


ANALYSIS OF THE DATA 

The data from these four satellites must be combined to form one overall picture of the artificial 
radiation belt. To do this assume that the energy spectrum of the electrons being counted is a fission 
spectrum. This is certainly the best guess. We will compare the data on this basis and see if there 
is agreement in the regions where direct comparison is possible. The fission energy spectrum N(E) 
is shown in Figure 1, curve A. A calibration of the Telstar I detectors at the Los Alamos Scientific 
Laboratory in a fission electron beam gives f, the fraction of fission electrons counted by the de- 
tectors, equal to 1/2.8 for the 240-340 kev channel and 1/6.0 for the 440-680 kev channel. 

For Injun we have the experimentally determined factor 1/f of several thousand, by comparison 
of two detectors on board. The 213 GM counter has also been calibrated at Los Alamos with a fission 
electron spectrum (private communication from A. Petschek, H. Motz, and R. Taschek, Los Alamos 
Scientific Laboratory), and the factor f determined this way is 1/4000. We will use this factor in the 
present analysis. The Los Alamos tests show that the detector counts bremsstrahlung from electrons 
of several Mev rather than direct penetrating electrons. (If the shield had been carbon rather than 
lead, the counter would have counted direct penetrating electrons.) 

For Traac, f is determined by considering the penetration of electrons through the detector shield 
of 0.265 gm/cm 2 of Mg and through the wall of 0.400 gm/cm 2 of stainless steel. Using the range 
straggling data (Reference 2) for A1 we can get the fraction of electrons that penetrate a shield of 
given thickness, as shown in Figure 2. The expression for the extrapolated range R is 

R = 0 . 526 E - 0.094 . 
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FRACTION (f) OF FISSION ELECTRONS COUNTED 



ENERGY (Mev) 

Figure 1— Curve A is the fission energy spectrum and curve B the transmission energy 
spectrum for the Traac GM counter (0.66 gm/cm 2 wall). 



ENERGY (Mev) 

Figure 2— The fraction of electrons of different energies that penetrate 
different shield thicknesses of Al. 
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This yields the absorber thickness that gives 10 percent transmission for electrons of energy E. For 
50 percent transmission we multiply the energy by 1.38, and for 80 percent transmission increase 
the energy by a factor of 1.92. In this way we get the electron transmission spectrum, curve B in 
Figure 1. The energies of the transmitted electrons are different from curve B, but the number 
transmitted is given correctly. The integral under this curve gives f = 1/5.5 for the Traac counter. 
More information on shielding calculations is given in Appendix A. 

Using the factors for the several detectors, we can calculate the total flux of fission electrons. 

In order to compare the different detectors, the total flux along several field lines (actually narrow 
ranges of L) has been plotted for different values of B (Figure 3). These plots show that the different 
detectors agree fairly well in flux values. Avoiding the first day after the nuclear explosion (labeled 
by the number 0 inside the symbols on the graphs) we can see quite smooth trends in the data. The 
flux from Telstar I may be as much as twice as high as Injun fluxes. Traac and Injun agree quite 
well where comparisons are possible. In general, the data shows agreement to a factor of 2. 

This agreement of the data shows two things: First, because the detectors give internally con- 
sistent results it seems likely that all the detectors are giving accurate information. Secondly, the 
assumption that the electrons have a fission energy spectrum appears to be correct. Of course it is 
possible that the energy spectrum is not a fission spectrum and also that the detectors are not in 
agreement, but it would have to be a peculiar combination of such effects that would give the agree- 
ment shown here. A comparison of the four channels of the Telstar I electron detector also indicates 
that the energy spectrum is fission-like up past 1 Mev. 


FLUX PLOTS 

Now that it has been demonstrated that the energy spectrum is essentially a fission spectrum at 
least in the region of data overlap we can use all the counter data to construct a composite flux map 
in B-L space. As Mcllwain has shown, these magnetic coordinates are the best way of organizing 
data about trapped particles (Reference 3). L is constant along a field line in space and, for a dipole, 
is the distance from the center of the earth to the equatorial crossing of the line, in units of earth 
radii (Figure 4). Values of L are calculated from the real values of the earth's field. 

In constructing the flux map for B > 0.15 gauss and for L < 2.0 earth radii the graphs in Fig- 
ure 3 are used to locate the flux contours. The experimental data outside this B-L region are essen- 
tially all from Telstar I. There are several weeks data from Telstar I and considerable redundancy. 
The map made this way is quite complete. The data available in early September gave the flux map 
in Figure 5. This map is for about 1 week after the explosion. There was considerably more flux at 
low altitudes at early times. 

This same data plotted in R-A coordinates, where 

M ,/ 3R 
B R3 r " L • 

R = L cos 2 A 
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B (gauss) 

Figure 3— The electron flux distributions along different field lines. 
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Figure 3 (continued)— The electron flux distributions along different field I 














B (gauss) 



gives an equivalent dipole representation of the 
earth's field (Figure 6). The maximum electron flux 
is about 2 x 10 9 elec/cm 2 -sec. Integrating to get the 
total number of electrons stored in the field we find 


\<p dV 


2 x 10 26 electrons. 


About 40 percent of these electrons lie inside the 10 9 
contour and about 60 percent of these electrons lie in- 
side the 3 x 10 8 contour. It is not certain what fraction 
of these electrons are bomb-induced and what fraction 


Figure 4— The B-L magnetic coordinate system. 

are natural electrons. In this region around L < 1.5 the energy spectrum seems softer than a fis- 
sion spectrum. 


The B-L flux map when plotted in terms of geographic coordinates gives the flux contours 
for different altitudes shown in Figure 7. 



Figure 5— The B-L map of electron fluxes. 
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AT 200 km 

Figure 7— Electron flux maps at different altitudes above the earth’s surface. Flux is in units of 10 5 

electrons/cm 2 -sec. 
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AT 400 km 

Figure 7 (continued)— Electron flux maps at different altitudes above the earth's surface. Flux is in units of 10 5 

electrons/cm 2 -sec. 
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AT 800 km 

Figure 7 (continued)— Electron flux maps at different altitudes above the earth's surface. Flux is in units of 10 5 

electrons/cm 2 -sec. 
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AT 1000 km 


Figure 7 (continued)-Eleclron flux maps at different altitudes above the earth's surface. Flux is in units of 10 

electrons/cm 2 -sec. 


VEHICLE-ENCOUNTERED FLUXES 

A machine code has been developed which calculates the total number of electrons/cm 2 in the 
artificial radiation belt that strike a vehicle in space. This is done by calculating a point on the vehi- 
cle trajectory, transforming to B-L coordinates, looking up the electron flux, and integrating along 
the vehicle orbit. This has been performed and the encountered fluxes have been determined for all 
of the vehicles listed in Table 1. These fluxes have been transformed into r/day by using 3 x 10 7 
electrons/cm 2 = lr. The orbital elements of these vehicle trajectories are given in Table 2. 

From the encountered fluxes in Table 1, we can learn several things. Let us first consider solar 
cell damage. The Bell Telephone Laboratories staff have studied this problem in considerable detail 
and prepared Figure 8, which shows how different type cells are damaged by 1 Mev electron irradia- 
tion (Reference 4 and private communication from the authors of that paper). Above about 0.5 Mev 
the electron damage is essentially independent of energy. Some care must be exercised in using this 
chart because of the variation in the characteristics of solar cells. We will assume all the electrons 
in the flux spectrum in Figure 5 are greater than 0.5 Mev in estimating the solar cell damage. 

About 20 percent degradation was needed by the blue sensitive p-on-n type cells on Ariel I to 
produce the observed power supply damage (private communication from A. Franta, Goddard Space 
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Table 1 

Calculations on Fluxes Encountered by Satellites Moving Through the Artificial Radiation Belt. 
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To be launched. 
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*To be launched. 

tSimilar to Nimbus, the Fixed Frequency Topside Sounder, the Swept Frequency Topside Sounder, and the Polar Ionosphere Beacon (which 
have been or will be launched). 
tSix orbits only. 




Table 2 


Orbital Elements of Various Space Vehicles. 


Element 

Ariel I 

Traac and 
Transit IV-B 

Telstar I 

Tiros V 

OSO I 

Relay 

Epoch 

(days, hours, 
min, sec) 

190, 9, 
0, 0 

190, 4, 
3, 46.506 

191, 8, 
51, 0 

190, 9, 
0 , 0 

190, 9, 
0 , 0 

305, 0, 
0 , 0 

Semimajor Axis 
(earth radii) 

1.1254 

1.1618 

1.5182 

1.1224 

1.0900 

1.5407 

Eccentricity 

0.05714 

0.009922 

0.2430 

0.02663 

0.003012 

0.2143 

Inclination 

(degrees) 

53.866 

32.423 

44.803 

58.102 

32.855 

50.0003 

Right Ascension 
of Ascending 
Node (degrees) 

-24.881 

96.434 

-156.222 

-75.536 

154.502 

163.708 

Argument of 
Perigee (degrees) 

-9.2537 

-51.6890 

164.811 

118.014 

139.136 

-167.526 

Mean Anomaly 
(degrees) 

-86.8833 

0.0001 

1.1684 

-194.11968 

-164.5453 

7.8219 


Table 2 (continued) 


Element 

Pogo 

1000 km 
Polar Orbit 

MA-7 

800 km 
Polar Orbit 

OAO 

SERB 

Epoch 

(days, hours, 
min, sec) 

82, 3, 
55, 32.101 

190, 9, 
0 , 0 

268, 14, 
0 , 0 

190, 9, 
0 , 0 

153, 0, 
0 , 0 

303, 10, 
0 , 0 

Semimajor Axis 
(earth radii) 

1.0931 

1.1568 

1.0331 

1.1254 

1.1270 

2.3284 

Eccentricity 

0.04830 

0.1490 x 10' 7 

0.008552 

4.4703 x 10' 7 

0.001074 

0.5518 

Inclination 

(degrees) 

90.001 

90.000 

32.546 

90 

30.982 

17.0 

Right Ascension of 
Ascending Node 
(degrees) 

-73.806 

-158.175 

75.069 

-158.175 

38.6643 

17.2908 

Argument of 
Perigee (degrees) 

-19.408 

-180.000 

78.188 

-180.000 

-68.8243 

134.6735 

Mean Anomaly 
(degrees) 

2.1956 

0.0000 

7.6908 

0 

36.7572 

° 
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Flight Center). This would be caused by about 
10 13 electrons/cm 2 according to Figure 8. 

About seven days after the nuclear explosion, 
this flux would have been achieved (Table 1 
gives 2.8 x 10 12 electrons/cm 2 -day for Ariel I, 
of which half hit the face of the cells). The 
Ariel I power supply started malfunctioning in 
3-1/2 days. This is quite good agreement. 

Traac and Transit IV-B also had blue sen- 
sitive p-on-n solar cells, but it would take 
3 x 10 M electrons/cm 2 to cause malfunction, 
because the cells were lower feficiency cells 
(private communication from R. Fischell, Ap- 
plied Physics Laboratory). Table 1 gives 
4.5 x 10 12 electrons/cm 2 encountered per day. 
Half of these electrons hit the face of the cells. 
Traac stopped transmitting in 36 days and 
Transit IV-B in 24 days. Using 30 days as the 
average, we get a total encountered flux of 
0.7 x 10 14 electrons/ c m 2 , only in fair agree- 
ment with that required to produce damage. 

Telstar I used the much more damage re- 
sistant n-on-p cells, because it was to routinely 
fly through the inner radiation belt protons. 

Even with the artificial radiation belt, its power supply lifetime is expected to be considerably 
longer than 1 year. 

The Telstar I solar cells are degrading at a rate that would be produced by 6 x 10 elec- 
trons/cm 2 -day of 1 Mev hitting the bare cells (private comrhunication from W. Brown, Bell Telephone 
Laboratories). This corresponds to about 1.8 x 10 13 electrons/cm 2 -day incident on the outside of 
the 30 mil sapphire covers. Our calculations give 1/2 X 2 x 10 13 = 1 X 10 13 electrons/cm -day 
hitting the cells. The observed solar cell degradation on Telstar I should be somewhat more than that 
calculated from the artificial electron belt, because slow proton damage probably contributes some- 
what to the degradation (private communication from W. Brown, Bell Telephone Laboratories). 

We have neglected the enhanced early time effects here on all the exposed satellites. An appreci- 
able part of the encountered flux may have been encountered in the first few days. For the first week 
after the explosion the flux was higher than that given in Figures 5 and 6. 

Injun, Tiros V, and other satellites continue to function. Injun has a low duty cycle and Tiros V 
shows some solar cell degradation. Film badge dosimeter measurements have been made on several 
space flights. About 10 r/day was measured inside 1.5 gm/cm 2 of shielding. In order to compare 



Figure 8-Solar cell damage curves. 
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this radiation dose with the predictions in Table 1, correction must be made for shielding. To 
do this we perform a calculation like that done for the Traac GM counter to get f , the fraction 


of electrons that penetrate the wall. Values of f 
have been calculated for different thicknesses of 
shield by using the relationship R = 0.526E - 0.094 
and the associated rough-straggling trans- 
mission curves in Figure 3. Figure 9 shows a 
plot of 1/f as a function of shield thickness. This 
is really only true for A1 but for lack of better 
information we will use it for other materials 
too. For 1.5 gm/cm 2 we get f = 1/50 for normal 
incidence particles. To correct for a distribution 
of incidence angles we will say roughly that about 
half as many get through. Also, 2n steradians 
are covered by a much thicker shield so that the 
total factor f' = 1/200. This would mean that 
10 r/day x 200 = 2000 r/day were incident on the 
outside of the vehicle. This agrees to within a 
factor of 2 with the calculated vehicle-encountered 
flux. 

MANNED FLIGHT 



SHIELD THICKNESS 


Figure 9— The fraction of fission electrons that penetrate 
different shield thicknesses. 


For a Mercury capsule orbit with an apogee of 264 km the total flux encountered in six orbits 
would be 0.24 x 10 10 electrons/cm 2 outside the vehicle (Table 3). If the apogee is lowered by 30 km 
(to 234 km) the total flux for 6 orbits is reduced to 0.17 x 10 10 electrons/cm 2 . If the apogee is 
raised by 30 km (to 294 km) the total flux for 6 orbits is increased to 0.45 x 10 10 electrons/cm 2 . 


PARTICLE TIME HISTORIES 

One of the important problems to answer 
about the new belt is how long it will last. The 
currently intense regions will last a number of 
years, according to present indications. At low 
altitudes the fluxes have already decayed a lot. 
According to Ariel I and Traac data, outside the 
10 5 contour of the B-L plot in Figure 6 the 
fluxes decayed several orders of magnitude in 
a few days. 


Table 3 

Flux per Orbit for a Mercury Capsule at an 
Altitude of 264 km. 


Orbit 

Flux (electrons/cm 2 ) 

1 

5.0 x 10 6 

2 

2.1 x 10 7 

3 

4.8 x 10 7 

4 

2.9 X 10 8 

5 

6.4 x 10 s 

6 

1.4 X 10 9 


Injun has noted some decay at 1000 km (private communication from B. O’Brien, State University 
of Iowa). At L = 1.18 and B - 0.191 there is a decay factor of about 2 from +10 to +1000 hr. 
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For the same L and time interval for B = 0.206 there is a decay factor of 4 and for B = 0.214 

there is a decay factor of 60. Injun saw no marked change in flux as a result of a modest size mag- 
netic storm. 

The only decay process we understand well enough to calculate is coulomb scattering. Particle 
time histories have been calculated for coulomb scattering and characteristic times determined (Ref- 
erence 5). The time to reach a scattering equilibrium (which is also about the time for this equilib- 
rium to decay to 1/e intensity) for different L values is listed in Table 4. Welch, Kauffman, et al. 
(Reference 5) first calculated these for solar maximum atmospheric densities and now, assuming that 
the density is less by a factor of 10, we get the values in Table 4. 

The air densities are not well known and the calculated times may be wrong by a factor of 5 or 
more. The calculated variation with L, however, should be fairly good. The Injun data show that the 
calculated times are of the right order of magnitude. The times show that the high flux region should 
last even through the next solar maximum if coulomb scattering is the principle loss process. 


Table 4 

Time Until Scattering Equilibrium for Different Values of L. 


L 

Calculated t 
( days) 

Measured r 

1.20 

10 

~ 1 month 

1.25 

150 

_ 

1.30 

1500 



1.35 

~3000 



1.40 

~10,000 

- 
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Appendix A 


Shielding and Radiation Doses 


Included here for the sake of completeness are some crude calculations on shieldings and dosages. 

One consideration that is important in some shielding calculations is bremsstrahlung. The doses 
delivered by the x-rays made by bremsstrahlung will be larger than the direct electron doses for 
large shield thicknesses. 

The fraction of the energy of an electron that goes into bremsstrahlung may be calculated from:* 

E b r e m ZE^ 

~Y. " 1600 - 

ion 

where z is the atomic number of the material involved. For the fission energy spectrum the average 
energy is about 1 Mev; 


F 

b rem _ C A 1 Fe Pb 

E ion ~ 0-004 0.008 0.015 0.050 


The energy spectrum of the x rays will be something like that in Figure Al. There will be a very 


few x rays up to 8 Mev, but not many over 2 or 
3 Mev. The low energy x rays (below about 100 
kev) will be absorbed in the shielding. This will 
remove about half the total energy in the x rays. 
The resultant transmitted energy spectra will 
have a peak at about 1/2 Mev (Figure Al). The 
x rays transmitted through the shield will be 
quite penetrating. Their mean free path will be 
roughly 20 gms/cm 2 . This means two things. 
First, they will be hard to absorb, and therefore 
it will take a lot more shielding to absorb them. 
Second, because they are hard to absorb, they 
will not be counted efficiently by a particle 
counter and also will result in less radiation 
dose. 

We can now calculate crudely the counting 
efficiency of the Injun (1961 o2) 213 GM counter. 

'Fermi, E., "Nuclear Physics,” A course given by Enrico Fermi 
at the University of Chicago, Notes compiled by Orear, 
Rosenfeld, A. H., and Schlurer, R. A., University of Chicago 
Press, revised edition, 1950. 



ENERGY (Mev) 


Figure Al— A crude bremssfrahlung x-ray 
energy spectrum. 
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From Figure 9 of the body of this report we see that it would only count about 1/20,000 of the fission 
electrons directly. But we find that 0.05 of the energy is converted to bremsstrahlung, of which half 
is absorbed in the shield. The mean energy of these x rays will be about 1/2 Mev. A normal GM 
counter will detect these x rays with about 1 percent efficiency. This gives 

( 0 . 05 )( 0 . 01 ) = 1/2000 , 

for the fraction of the electrons counted via bremsstrahlung. This calculation is not very accurate 
but it does show that the Injun counter counts electrons via bremsstrahlung with about the observed 
efficiency. 


Manned Flight 

The effects of the new radiation belt on manned flights must be considered. For the Mercury 
project the total flux that would be encountered for a six orbit mission with the MA-7 (1962 rl) orbit is 
0.24 x 10 10 electrons/cm 2 outside the vehicle, or 80 r (3 x 10 7 electrons/cm 2 = lr). The shielding 
of the vehicle is such that about 1 percent of this dose is delivered to the astronaut, about lr. 

The mercury dose is almost all received in the South Atlantic "hot spot" (see Figure 8 of the 
body of the report) and occurs mainly on orbits 4, 5, and 6. The breakdown of the lr dose inside the 
capsule by orbits is given in Table Al. 


Table Al 

Radiation Dose per Orbit for a 
6 Orbit Flight on the MA-7 Orbit. 


Orbit 

Dose (r) 

1 

0.003 

2 

0.01 

3 

0.03 

4 

0.1 

5 

0.3 

6 

0.6 
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